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Summary. CoUapsars may be a source for the "long" Gamma Ray Bursts (GRBs) 
in the BATSE catalog. CoUapsars may radiate gamma rays anisotropically by 
beamed jet emission close to the observer's line of sight. These jets must penetrate 
the radiation-dominated medium of their collapsar progenitor and break through its 
atmosphere in order to produce a GRB. We present a study of previously published, 
axisymmetric, collapsar jet simulations. Here we use the linearized, relativistic fluid 
equations to find the Kelvin-Helmholtz modes that are triggered by recollimation 
shocks within the jet. The modes will grow as they propagate with the jet. These 
are of interest since the light curves of GRBs show evidence of a variable flow from 
the GRB engine. We also evaluate effects of grid scaling in the numerical simulation 
and show that short wavelength modes are suppressed by grid scaling before the jet 
breaks out of the Helium shell. 



1 Introduction 

The close association between GRB030329 and SN2003dh proves that a subset 
of "long" GRBs is associated with core collapse supernovae [1,2]. The light 
curve of GRB030329 steepens over time. Price et al. [3] interpret this as a 
broadening of a jet-like outflow from the GRB source. Multiwavelength obser- 
vations by Uemura et al. [4] along with those from GRB021004 [5] strengthen 
the case for a variable jet outflow from the central engine. 

These observations support the prompt (Type I) collapsar model of GRBs 
[6,7]. In this model, a progenitor star with mass > 30 Mq fails to eject its 
envelope after the core collapses into a black hole. Accretion of free falling 
material may drive a relativistic jet of gas that may break through the stellar 
envelope. At issue is whether hydrodynamic instabilities within the jet will 
grow so that the radial profile of jet velocity becomes time and space depen- 
dent. Such information could offer insight into the production of gamma rays 
as well as prompt optical afterglows of GRBs. 
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2 Dispersion Relations for Collapsar Jets 

Analyses of collapsar jet models have shown that the jet is unstable as it prop- 
agates through its progenitors atmosphere. Axisymmetric, relativistic hydro- 
dynamic simulations by Aloy et al. (2000) [8] show evidence of pinch-body 
modes triggered by recoUimation shocks. 

We studied the simulations that Aloy et al. (2000) [8] made using the au- 
thor's GENESIS code. In the radial direction, the energy deposition region 
extends from the inner grid boundary located at 200 km to a radius of 600 
km (2-6 X 10^ cm) . We studied two of their jets with constant energy depo- 
sition rates of dE/dt = 10^° (C50) and 10^^ ergs s~^ (C51). Our approach to 
analyzing the space-dependent structures is to use the linearized relativistic 
fluid equations to find the Kelvin-Helmholtz modes that may be triggered 
by recoUimation shocks within the jet. Space-dependent perturbations prop- 
agate with a dispersion relation given in the appendix of Hardee, Clarke & 
Rosen (1997) [9]. This dispersion relation assumes uniform conditions inside 
and outside a jet with sharp discontinuity at the jet surface. 

The jets in the simulations have velocity profiles in the transverse direction. 
In order to identify modes in the simulation we must determine typical sound 
speeds in the jet and the external medium. To calculate the value of the 
relevant parameters in the jet and the external medium, we define the jet 
radius Rj wherever Fvr/c = 0.5r*v*/c. Here Vr is the radial velocity, F = 
Xj \J\ — vll (? is the Lorentz factor, v* is the maximum value for the jet 
velocity at a given cross section and F* the corresponding Lorentz value. 
Figure 1 shows the angular cross sections of pressure, density and velocity 
for the jet C50 near its base and near the head of the jet. We calculate the 
weighted average for jet speed and sound speed in the region R < Rj and 
sound speed in the interval Rj < R < 3Rj. The local sound speed for a 
relativistic gas is defined by 



where 7 is the adiabatic index (4/3 for a relativistic gas), p is the pressure, p 
is the density, and c is the speed of light. From the average of the jet speed 
typical values of F range from 1-5; however, the C51 jet late in the evolution 
near the jet head reaches a. F ^ 20. 



3 The Stages of Jet Evolution 




(1) 



We identify two distinct stages of jet evolution: downstream from the inlet and 
upstream from the breakout. The inlet region corresponds to where the jet is 
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model 50 3.35s r = 6.5 x 10' cm 


model 50 3.35s r = 1.3 x 10'° cm 
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Fig. 1. Cross sections for pressure P, density p, velocity v^ , and sound speed a^ for 
the C50 jet at t=3.35s for cuts near tlie base and tlie jet head. Black vertical lines 
identify the defined jet radius, Rj. 



injected. In the C50 jet simulation at time t=3.35 s, we can study the stability 
downstream from the inlet near its base (r ~ lO^cm). Breakout through the 
Helium shell of the progenitor occurs at r ~ lO^^cm. Here the stability can be 
studied in the region immediately upstream from the jet head (r ~ 10^*^0771). 
The dynamics of these two regions are comparable to the same jet at a later 
time and for the C51 simulation. Figure 2 shows the result of our stabil- 
ity analysis for these two regions with our calculated values for normalized 
wavelength and frequency for the fundamental, 1st, 2nd, and 3rd pinch body 
modes. 



Once the jet has penetrated the He shell, it enters a hydrogen medium 
where it will be overpressured. The jet will expand adiabaticaly in the per- 
pendicular direction to the jet axis at cj^fZ in the reference frame of the 
fluid [10,11]. RecoUimation shocks below the He-H boundary shell destroy 
all upstream perturbations. At the same time they can trigger pinch modes 
downstream. We are most interested in perturbations just upstream of the jet 
head after it has crossed the He layer at lO^'^cm. It is at this point that the 
jet becomes a causally disconnected wind where density and velocity pertur- 
bations "freeze" in the base of the wind. 
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Fig. 2. Top row diagram shows cuts for the C50 jet at t=3.35 seconds for the jet 
close to the base (6.5 x 10* cm) and for the jet close to the head (1.3 x 10^" cm). 
Cuts are for the angles 0.25°, 1.125°, 2.25° and 3.25° (solid, long dash, .short da.sh 
and dash-dot correspondingly) . On the bottom row diagrams show the fundamental 
(SO), 1st, 2nd, and 3rd pinch body modes (Bl, B2, B3). These are the real and 
imaginary solutions for these regions. 



4 Conclusion 



We solved the dispersion relations in Hardee, Clarke, & Rosen 1997 [9] for 
fc(w) with the appropriate parameters from each cut of the jets. We found 
the solutions for the fundamental and the first three body modes for both 
jets. Figure 2 shows the solutions for C50 just downstream from the inlet 
and upstream from breakout at t=3.35s. We compute the wavelengths at 
the maximum growth rate A* = 2t:vI^ /uj* (where and w* are the wave 
speed and frequency at maximum growth) and the minimum growth length 
I* = {k^'^ (where k*j is the imaginary wave number). The growth time 
T = i/vw for any given solution is calculated from. 



k\ + k'j 



(2) 
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where kfi is the real wave number. We find that our solutions depend on 
the distance from the inlet r. The range for maximum growth wavelength A* 
for the body modes near the jet inlet is O.lr to r, and the range upstream 
from breakout is 0.06r to 0.6r. The range for the growth length £* is 0.5r to 
r near the inlet, and for near breakout the range is 0.2r to 2r. Near breakout 
the growth time r ranges from 0.3 r/c to 3 r/c or 1 to 10 seconds at the 
distance of the He shell. 

Our numerical analysis of jet stability assumes a sharp lateral boundary 
between the jet and the external medium. The analytical approach of Aloy et 
al. (2002) [12] for the same jet asumes an extended shear layer. For Lorentz 
factors of order 10, they predict a growth time r ~ O.Olr/c which is one order 
of magnitude smaller than our method; however, the characteristic order of 
A* ^ O.lr — 0.5r is also what we predict. 

One objective of our study was to find if there were wave modes that would 
not show up in the simulation because their length scale would be smaller than 
the computational grid scale. Aloy et al. (2000) use a logarithmic grid scaling 
so that the jet is better resolved at the base than at the head of the jet [8]. 
For the jet downstream from the inlet, grid scaling is ^ 0.1 Rj and this is 
sufficient to resolve wave modes. Near the jet head the grid scaling is ~ Rj 
and this can suppress the pinch body modes numerically. Since this region will 
be the most relevant to the evolution of the jet flow at breakout, simulations 
with greater grid resolution need to be done. 

E. A. G. wishes to acknowledge the funding of this study through the 
Alabama Space Grant Consortium project. The authors would like to thank 
Miguel Angel Aloy for his generosity in sharing his simulation results and his 
comments. 
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